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We further elaborate the theory of quantum fluctuations in vertical-cavity surface- 
emitting lasers (VCSELs), developed in Ref. In particular, we introduce the 
quantum Stokes parameters to describe the quantum self- and cross-correlations 
between two polarization components of the electromagnetic field generated by this 
type of lasers. We calculate analytically the fluctuation spectra of these parameters 
and discuss experiments in which they can be measured. We demonstrate that in 
certain situations VCSELs can exhibit polarization squeezing over some range of 
spectral frequencies. This polarization squeezing has its origin in sub-Poissonian 
pumping statistics of the active laser medium. 
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I. INTRODUCTION 

In the last years there has been an increasing interest to the polarization properties of 
the vertical-cavity surface-emitting lasers (VCSELs). This interest is motivated in the first 
line by the potential applications of this type of lasers in the high-rate optical communica- 
tions But there is also more fundamental reason for understanding of the polarization 
behavior in VCSELs, namely, a possibility of generating the intensity-squeezed light using 
the sub-Poissonian pumping of the active medium 0,0]. To date, squeezing in VCSELs has 
been demonstrated experimentally for both single-mode operation and in a multi-transverse- 
mode regime In single- mode operation with only one linearly polarized mode above 
threshold the fluctuations in a sub-threshold mode with polarization orthogonal to the lasing 
mode can present large intensity noise 0, Q and, moreover, be highly correlated with the 
intensity fluctuations of the oscillating mode. This phenomenon can result in deterioration 
of squeezing observed in experiments with polarization sensitive optical elements. Therefore, 
polarization dynamics in VCSELs plays an important role for correct description of their 
quantum fluctuations. 

At present, the standard theory that accounts for dynamics of two polarization compo- 
nents of the electromagnetic field in VCSELs is the so-called "spin-flip" model developed 
by San Miguel, Feng and Moloney j^. On the basis of this model several authors have 
formulated semiclassical theories of light fluctuations in VCSELs f3[ 0, [l^ l], [3| • How- 
ever, semiclassical description is inappropriate for intensity-squeezed light and, therefore, 
calls for fully quantum model of quantum fluctuations in VCSELs. The "quantum spin-flip" 
model was developed recently in Ref. This model takes into account, on the one side, 
the dynamics of two polarization components of the electromagnetic field and, on the other 
side, the pumping statistics of the active laser medium. In particular, the quantum spin-flip 
theory allows for sub-Poissonian pumping statistics in which case VCSELs can generate the 
intensity-squeezed hght. 

In this paper we further elaborate the quantum spin-flip model of VCSELs, developed 
in In particular, we apply the quantum Stokes parameters to describe the quantum self- 
and cross-correlations of two polarization components of the electromagnetic field generated 
by VCSELs. We analytically calculate the fluctuation spectra of the quantum Stokes pa- 
rameters and discuss experiments in which they can be measured. We demonstrate that for 
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the sub-Poissonian pumping statistics VCSELs can exhibit polarization squeezing in some 
range of spectral frequencies. 

The paper is organized as follows. In Sec.|n]we give a short resume of the quantum spin- 
flip model developed in Ref . Q] , and calculate analytically the spectral densities of quantum 
fluctuations of the quadrature components. In Sec. IIIII we introduce the quantum Stokes 
parameters, their fluctuation spectra, and discuss the experiments where these fluctuations 
spectra can be measured. Using the results obtained in Sec. HT] we analytically calculate 
the fluctuation spectra of the quantum Stokes parameters. In Sec. IIVI with help of the 
analytical results obtained in Sec. lIIII we illustrate graphically the possibilities of observation 
of polarization squeezing in VCSELs. We also provide the figures of typical cross-correlations 
spectra of photocurrents and cross-correlation spectra of the Stokes parameters 5*2 and 5*3 
that can be measured experimentally. In Sec.|3we summarize the results. 



II. QUANTUM SPIN-FLIP THEORY OF VCSELS 



A. Resume of the model 



In this section we shall give a brief resume of the quantum spin-flip model of VCSELs 
developed in Ref. jl| . We shall define the physical parameters of this model and provide the 
equations which will be used in the following sections. For more details we refer the reader 
to Ref. Q. 

The semiclassical four-level spin-flip model of VCSELs was developed by San Miguel, 
Feng and Moloney j^. This model describes very well the dynamics of these semiconductor 
lasers and is widely used for understanding of such phenomena, for example, as polarization 
switching. The spin-flip model takes into account the spin sublevels of the total angular mo- 
mentum of the heavy holes in the valence band and of the electrons in the conduction band. 
These four sublevels interact with two circularly polarized electromagnetic waves in the 
laser resonator and it is this interaction that is responsible for the complicated polarization 
dynamics manifested by this type of lasers. 

The four- level scheme of the semiconductor medium is shown in Fig. L Two lower levels 
\b±) correspond to the unexcited state of the semiconductor medium with zero electron-hole 
pairs while the upper levels |a±) to the excited states with an electron-hole pair created 3|- 
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Two pairs of levels |a+), \b+) and \a—), \b—) are coupled via interaction with the left and 
right circularly polarized electromagnetic waves in the laser cavity described by the field 
operators d+{t) and a_(t). As explained in Ref. j^, physically these two pairs of transitions 
are associated with two 2;-components Jz = ±1/2 of the total angular momentum J = 1/2 
of the electrons in the conduction band and corresponding z-components Jz = ±3/2 for 
J = 3/2 of the heavy holes in the valence band. The constants 7^ and 7^ are the decay 
rates of the populations of the upper and lower levels, 7^ (not shown in Fig. 1) is the decay 
rate of the polarization, and 7c is the spin-flip rate that accounts for mixing of populations 
with opposite values of Jz- The last parameter was introduced in Ref. to describe the 
spin-flip relaxation process. This parameter is responsible for coupling of two transitions 
with different circular polarizations and, as a result, for various polarization dynamics of 
VCSELs. 

It should be noted that the authors of Ref. have considered the situations of equal 
relaxation constants of the upper and the lower levels, 7a = 7b. However, it is known 

nn 

from the literature y, m that this is not the most favorable condition for generation of the 

n 

sub-Poissonian light. Therefore, the quantum spin-flip theory in Ref. y was developed for 
arbitrary values of 7^ and 7;,. In this paper we shall also consider this general situation. 

Moreover, it has been mentioned in the literature (see, for example, Ref. that this 
model describes correctly a semiconductor laser if we assume the decay rate 7^ of the lower 
levels to be very large compared to the other decay constants, namely, 7a, 7c and k. From 
the classical point of view both situations 7;, = 7a and 7^ ^ •ja result in the same dynamical 
behavior of VCSELs. However, it turns out that the statistical properties of two models 
with 7f, = 7a and % ^ 7a are very different. The detailed discussion of this difference is out 
of the scope of this paper and we shall address this point elsewhere. 

We have indicated in Fig. 1 the pump process with mean total pumping rate 2R which is 
then separated with equal probabilities between two sublevels |a+) and |a_). Quantum spin- 
flip model of Ref. jl| takes into account a possibility of sub-Poissonian pumping of the laser 
medium using the technique of the pump- noise suppression 0, 0]. For stationary in time 
average pumping rate, the influence of the pump statistics can be characterized by a single 
parameter p < 1 jl^. Iisj]. For p = 1 the pump is perfectly regular while for p = the pump 
has Poissonian statistics. Intermediate values of < p < 1 correspond to sub-Poissonian 
pumping while for p < the pump process possess the excess classical fluctuations and 
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corresponds to super-Poissonian statistics. 

This pump statistics was introduced into the quantum spin-flip model using the 
Heisenberg-Langevin equations for the operator- valued collective populations Na±{t), Nb±{t) 
of the upper and lower levels in Fig. 1, and for the collective polarization P±(t). On the 
basis of the Heisenberg-Langevin equations, the equivalent c-number Langevin equations 
were derived for the collective atomic and field variables, corresponding to the normal or- 



dering of the atomic and field operators 



15l |. Next, using the fact that the relaxation 



rates 7b of the lower levels and 7^ of the polarization in VCSELs are much bigger than the 
relaxation rate 7a of the upper levels, the macroscopic c-number populations Nb±(t) and 
the macroscopic c-number polarization P±{t) were adiabatically eliminated. The resulting 
equations can be written in terms of the total population of two upper levels |a+) and 
|a— ), and of the total inversion between them. The corresponding variables are defined 
as D{t) = (Na+it) + Na-{t))/2, and d{t) = (Na+it) - Na-{t))/2. The equations for these 
variables and the two c-number field components a±{t) are 

d±(t) = -Ka±{t) - + iuJp)a^{t) + c(l - ia)\p{t) ± d{t)\a±{t) + (2.1) 

b(t) =R- ^D{t) - c(|a+(t)|2 + \a_{t)\^)D{t) - c(|a+(t)p - |a_(t)nrf(t) + (2.2) 

d{t) = -^,d{t) - c(|a+(t)p - |a_(t)nD(t) - c(|a+(t)|2 + \a^{t)\^)d{t) + Fait). (2.3) 

Here /t is the cavity damping constant, Up and describe the linear birefringence and the 
linear dichroism of the semiconductor medium. The last parameter was not included into 
the model in Ref. |Q and is introduced here as a generalization. Next, a is the linewidth 
enhancement in semiconductor lasers, 

a = — , (2.4) 

where v is the frequency of the semiconductor energy gap, and uj is the resonator frequency. 
We have also defined the relaxation rate 7^ as 7s = 7a + 27c, and have introduced the 
following shorthands, 

7^(1 + a^) 

where g is the coupling constant of interaction of the electromagnetic field with the polar- 
ization. 

The functions F±{t), Fu{t) and Fd{t) are the c-number Langevin forces. Their nonzero 
correlation functions were calculated in Ref. Q . In general the results are rather cumbersome 
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but they are simplified in the vicinity of the stationary solutions. For completeness we shall 
give these correlation functions for the stationary solutions at the end of this section. 

B. Stationary semiclassical solutions 

Semiclassical equations of VCSELs are obtained from Eqs. ()2.1|) - ()2.3p by dropping the 
c-number Langevin forces. In this subsection we shall give the stationary solutions of these 
equations which characterize the stationary generation of VCSELs. For investigation of 
quantum fluctuations in VCSELs we shall use standard assumption that these fluctuations 
are small compared to the corresponding stationary values. This will allow for linearization 
of Eqs. ()2.H) - ()2.3|) around stationary solutions with respect to the quantum fluctuations. 

Stationary solutions of Eqs. ()2.H) - ()2.3p have been investigated in detail in 0, [l^. When 
Up and 7^ there are in general four types of stationary solutions: two of them have 
linear polarization along the x and y axes, and two other elliptical polarization. We shall 
consider only linearly polarized solutions because this type of solutions is usually realized in 
experiments. In this case the circularly polarized field components have equal amplitudes 
and can be written in the form 

a±it) = Qe'^^'^^\ (2.6) 

where the real amplitude Q is normalized so that = | a+ P = | ct- P gives the mean number 
of photons in the corresponding circularly polarized field mode. Two other parameters A 
and determine the type of polarization of the stationary solution ()2.6|) . 

We remind that the linearly polarized field components a^it) and ay{t) are related to the 
circularly polarized ones as 

,,,(i) = M«H^, = (2.7) 

For the x-polarized solution = and for the y-polarized solution = 7t/2. The frequency 
detunings A in Eq. ()2.6|1 are different for these solutions and are equal to 

'^x,y = - [l^x.yOt ± UJp] , (2.8) 

where the upper sign corresponds to the x-polarized solution and the lower sign to the y- 
polarized one. Here we have introduced the shorthands = + t^a and Hy = k — Ha- The 
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x-polarized stationary solution reads 

= \f2Qe'^^\ ay = 0, (2.9) 
while the ^/-polarized stationary solution is given by 

= 0, ay = V2Qe'^yK (2.10) 

For both solutions we have 



Q = ^/,(r-l), (2.11) 

where r = R/Rth is the dimensionless pumping rate, i?th is the threshold pumping rate, and 
Is is the saturation intensity; the two latter are given by 

i?th = ^, Is = ^. (2.12) 
c 2c 

Note that for > the threshold pumping rate for the y-polarized solution is lower that 
for the x-polarized one. 

The stationary values of the atomic variables do and Dq for these linearly polarized 
solutions are equal to 

^0 = 0, ^0 = —^ = ^. (2.13) 
7 + 2cQ^ c 

In the case of VCSELs as in general for solid-state and semiconductor lasers the question of 
stability of stationary solutions is very important. The stability analysis of these stationary 
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17l |. and 



solutions was performed in a number of publications, as for example, Refs. 
we refer the reader to this papers for details. In our analysis of quantum fluctuations we 
shall assume that the corresponding stationary operation regime of VCSEL is stable. Since 
for low pumping rate only x-polarized solution is stable, we shall restrict our analysis of 
quantum fluctuations only for this type of stationary solutions. 



C. Linearization around stationary solutions 

To calculate the quantum fluctuations around the stationary solution we shall linearize 
Eqs. ()2.H) - ()2.3|) around the steady state given by Eq. ()2.6|) . As mentioned above we shall con- 
sider here only x-polarized stationary solution. Adding small fluctuations to the stationary 
solutions we can write the field and the atomic variables as 



a±{t) = {Q + Sa±{t))e'^\ D{t) = Dq + 6D{t), d{t) = 6d{t). (2.14) 
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In this equation and in what follows we have dropped the index x in A^, since we shall be 
concerned only with x-polarized solution. Substituting these expressions into Eqs. fj2.1|) - ()2.3|) 
and linearizing, we arrive at the following equations for small fluctuations, 

^5a±{t) = {Ka + iup)(5a±{t) - da^it)) + c{l - ia)Q{5D{t) ± 5d{t)) + F±{t)e-'^\ 
j^5D{t) = - (7 + 2cg2) 5D{t) - K.,Q (Sa+it) + 5a_(t) + c.c.) + ^^(t), 
j^5d{t) = - (7s + 2cQ^) Sd{t) - K^Q {5a+{t) - 5a_{t) + c.c.) + Fd{t). (2.15) 

It is convenient to introduce the fluctuations of the linearly polarized components of the 
field Sa^(t) and 6ay(t), defined according to Eq. ()2.7|) . for which the set of coupled equations 
()2.15|) decouples in two sets of independent equations for Sa^it) and 6ay(t) with Langevin 
forces F;j.{t) and Fy{t) defined similar to Eq. ()2.7|) . Moreover, we shall define the fluctuations 
of the amplitude and the phase quadrature components, SX^it) and SYx{t) of the a;-polarized 
field component, 

6X,{t) = ^{6a^{t) + 6a:{t)), 5Y,{t) = ^((5a,(t) - KW), (2-16) 

and similar for the ?/-polarized component. For these fluctuations we obtain the following 
equations, 

j^SX^it) = V2cQ6D{t) + R,{t), 
j^^Yxit) = -V2acQ6D{t)+T,{t), 

j^5D{t) = -T6Dit) - 2V2K,Q6X^{t) + Foit), (2.17) 

and 

Y^Xyit) = 2KjXy{t) - 2tUp5Yy{t) - V2acQ5d{t) + Ry{t), 

(JjL 

f^SYyit) = 2KjYy{t) + 2iUpSXy{t) " V2cQ S d{t) + Tyit), 
jM{t) = -VMit) + 2V2K,Q5Yy{t) + Frf(t), (2.18) 

where the new Langevin forces Rx{t) and Sx{t) are defined as 

R.it) = i(F.(t)e-^* + F;(t)c^^*), Tx{t) = ^(i^.(t)e-^^* - F:(t)e^^*), 

Ry(t) = l{Fy{t)e-'^' + F;{t)e^^% Ty(t) = ^(F,(t)c-^^* - F;(t)e^^*). (2.19) 



9 



In Eqs. fj2.17|) and ()2.18|) we have introduced 

r = 7 + 2cQ2 = 7r, = 7, + 2cg2 = 7, + 7(r - 1), (2.20) 
as convenient shorthands. 

D. Spectral densities of quantum fluctuations 

To solve Eqs. ()2.17|) and ()2.18|) we take the Fourier transform of the field and atomic 
fluctuations, 

^ +00 

SX^{n) = ^= f 5X^{t)e^^^dt, (2.21) 

and similar for the other variables, that converts these differential equations into algebraic 
ones. The spectral correlation functions of these quadratures are (5-correlated, 

{5X,{n)5X,{n')) = {5X^)^5{n + n'), 
{5Y,{Q)5Yi{n')) = {5Y^)nm + ^'). 

{6X,{Q)5Y,{n')) = {6X,6Yi)n6{n + n'), (2.22) 

with (5Xj^)n, i = x,y and {SY^^)q being the spectral densities of the corresponding quadra- 
tures, and {6Xi6Yi)^ their cross-spectral density. 

After a simple algebra we obtain the following expressions for the fluctuations of the 
amplitude quadratures SX^i^) and 6Xy{Q), and the phase quadrature 6Yy{Q): 

sx,{n) = ;^^{(r - tn)R,{n) + V2cQFn{n)}, (2.23) 
sxy{n) = -^i^{[2/€,7(r - 1) - {2Ka + in){Ts - in)]Ry{n) 

- [2aK^-f{r-l) + 2ujp{r,-in)]Ty{Q) + V2cQ{2up + 2aKa + ian)Fdin)}, (2.24) 

^^'^^^ = -^{2^,(r,-^(])i?,(f]) 

- {2Ka + in){Ts-in)Ty{n) + V2cQ{-2aujp + 2Ka + in)Fd{n)}, (2.25) 

with 

D^{n) = -in{T -in) +2K,^-f{r -1), 

Dy{Q) = (r, - iQ)[{2ujpy + {2Ka + i^lf] + 2K^7(r - l){2aujp - 2Ka - iVt). (2.26) 
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The other phase quadrature SYx{Q) will not appear in the observables that we shall discuss 
below. Using the results obtained in Ref. |^and taking into account the stationary solu- 
tions ()2.(ij) and ()2.13|1 we obtain the following nonzero correlation functions of the Langevin 
forces Ri(t),Ti(t) with i = x,y, and FD(t),Fd(t) for the stationary regime of VCSEL in 
approximation of the small fluctuations, 



{R,{t)R,{t')) = {Ry{t)Ry{t')) = {T^{t)T^{t')) = {Ty{t)Ty{t')) = K J {t - t' ) , 

{Fn{t)Fn{t')) = ^t{i ~\p)6{t - t% 

{Fa{t)F,{t')) = ^T,5{t-t'), 
c 

{FD{t)R^{t')) = {Fait)Ty{t')) = -V2K,Q6{t-t'). 



(2.27) 



Equations (j2.23|) - (j2.26p together with correlation functions (|2.27|) allow us to evaluate an 
arbitrary correlation function of the laser light emitted by the VCSEL. The spectral densities 
of the amplitude quadratures {6X^)^, {6Xy)n are given by, 



ml) 



2\Dy{n)\'' 

with Ax and Bx determined as. 



{fi2 + ^V[l-(r-l)p/2]}, 



[O.^ + AxVt' + ABx}, 



(2.28) 
(2.29) 



A 



X 



+ 



2ui!p + a7(r — 1) — 4/t7(r — 1 



B 



X 



"I 2 r "I z (-) 

Ka'^s- K.'lir -I) + idp'^s + oiK'^{r - I) + 7s7(r - 1) (ttK^ + c^p) , (2.30) 



The spectral density of the phase quadrature component (5K^)n is equal to. 



2\Dym^ 



f2.311 



with Ay and By given by. 

Ay = Ainl + ul) + 7^ + 7(r - l)(4a^p + 7s), 



B 



Y 



YsK + ^p) + Islir - 1) i^'pia' + 2) + < + u'pYir - lYia' + 1), (2.32) 



Finally the cross-spectral density {6Xy6Yy)fi reads, 
-Kxl{r - 1) 



(SXySYy), 



2\Dym' 



-{aKxn^+2Kujp-f{r~l){a^+l)+2-f, 



K{aKa + UJp) + aKa{Ka-aUJp) |, 

(2.33) 
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These analytical results will be used below for evaluation of the spectral densities of the 
quantum Stokes parameters, their cross-spectral densities and for the cross-correlation spec- 
tra of the photocurrents. 



III. QUANTUM POLARIZATION STATES OF LIGHT: GENERAL 

DISCUSSION 



A. Quantum Stokes parameters 



There are two equivalent descriptions of the polarization properties of light in classical 
optics either by the polarization matrix or in terms of the classical Stokes parameters jisl ]. 
During the last decade the quantum-mechanical version of the classical Stokes parameters 
was introduced in the literature and very actively used in quantum optics to describe the 
quantum fluctuations of polarization of the electromagnetic field 19, 20, 21, 2^. There have 
been several theoretical proposals for generation of polarization-squeezed li ght j2lL l23l . I24L 

0, 13 ^ experiments in which such kind of light was observed boi Isfli l8l|. 

We shall use the language of the quantum Stokes parameters for characterization of the 
quantum fluctuations of polarized light in VCSELs. In this section we shall express the 
fluctuation spectra of the quantum Stokes parameters through the spectral densities of the 
quadrature components evaluated above. In the next section we shall apply these results for 
the particular case of VCSELs. 

Let us write the operator E{t) of the electromagnetic field at the output of the VCSEL 
in terms of the x- and y-polarized components. 



E{t) = a^{t)e^ + ay(t)ey, 



(3.1) 



where a^it) and ay{t) are the photon annihilation operators in the Heisenberg representation. 
In what follows we shall omit the time argument when this does not create ambiguities. 
The quantum Stokes operators Sa,fi = 0,1,2,3 are introduced similarly to their classical 
counterparts (see, for example j2^). 



^0 

^1 

^2 
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Ss = i{aldx - aldy) . (3.2) 
Using the commutation relations for the photon annihilation and creation operators, 

[di,d]] = 6ij, {ij = x,y), (3.3) 
it is easy to verify that the operator 5*0 commutes with all the others, 

[So,S^]=0, (/i = 1,2,3), (3.4) 

and that the operators 5*1, 5*2 and 5*3 satisfy the commutation relations similar to the com- 
ponents of the angular-momentum operator, 

[S,, S2] = 22^3, [^2, 53] = 2tSi, [S3, Si] = 2iS2. (3.5) 

The noncommutativity of these three Stokes operators does not allow their simultaneous 
measurement in any real physical experiment. The mean values (S'^),/i = 1,2,3 and the 



variances AS*^ = y {{S^ — {S^)y) are given by the uncertainty relations 

A5iA52 > 1(4)1, AS2A53 > |(5i)|, 5S,5Sr > \{S2)\. (3.6) 

When the x- and ?/-polarized components of the electromagnetic field are in coherent states 
\ax) and \ay) i. e., 

one can speak about the coherent polarization state of the electromagnetic field. The mean 
values of the quantum Stokes parameters in this state are obtained by replacing d^ — > 
and dy ay in Eq. ()3.2|) . For example, for the first two parameters one obtains, 

(5*0) = + \ay\^ = {h^) + {fly) = (h), 

{Si) = \a^\^ - \ay\^ = (n^) - {fiy), (3.8) 

where (n) is the mean total number of photons in the electromagnetic wave. The variances 
of all four quantum Stokes parameters in this case are equal and given bv |27|. 

ASl = {n,) + {ny) = {n), /i = 0,l,2,3. (3.9) 

This property of the coherent polarization state allows one to define a polarization squeezed 



state similar to the definition of a single-mode squeezed state. According to 



2l| one can 
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speak about polarization squeezing if one of the four variances AS*^ of the Stokes parameters 
becomes smaller than that in the coherent state, i. e. AS*^ < (h) for at least one /x. 

Classical Stokes parameters Sfj_,fi = 0,1,2,3 (without hats) are obtained as the mean 
values of their quantum versions defined in Eq. ()3.2j) . S*^ = (S^). From the classical point of 
view, all polarization properties of light are completely described by these four parameters: 
5*0 determines the total beam intensity, while three other parameters characterize the polar- 
ization state of the light beam. This polarization state in classical optics is often represented 
in a Poincare sphere with Si, S2 and 5*3 forming its three orthogonal axes. 

In quantum optics to completely characterize polarization properties of light in addition 
to the mean values 5*^ of the quantum Stokes parameters one has to determine their variances 
AS"^. In general all these variances can be different and one can speak of an uncertainty 



ellipsoid in the Stokes-Poincare space 



22j | . In general case, when different Stokes components 



are correlated, there are three additional parameters which determine the orientation axes 
of this uncertainty ellipsoid. 

While the general description is outside of the scope of our paper, we shall illustrate below 
graphically that in the case of VCSELs different quantum Stokes components S*^ can have 
different variances AS*^. The quantum fluctuations of polarization in VCSELs are therefore 
characterized by an uncertainty ellipsoid in the Stokes-Poincare space. 



B. Measurement of the classical Stokes parameters 

Four classical Stokes parameters 5*^ can be measured in an experimental setup shown 
in Fig. 2. This measurement scheme consists of a compensator, a polarizing beam splitter 
(PBS), and two photodetectors. Let 6x and 6y denote the phase changes produced by the 
compensator in the x- and y-components of the electromagnetic field given by Eq. 1)3.11] . 
Next, let (fi denotes the angle between the transmission axis of the PBS and the x-axis. 
Then the field amplitudes di and 0,2 of the transmitted and reflected waves after the PBS 
can be written as 

CLi = e*^'' (d^. cos (f + dyC"'^ sin (p) , 

0.2 = e*''"' (— da; sin -|- dj^e~*^ cos 9?), (3.10) 
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where 6 = 5^ — Sy is the phase difference between the x- and y-components introduced by 
the compensator. 

The secondary waves after PBS are photodetected and one observes the mean values of 
the photocurrents {ii) = 7]c{a\ai), and {12) = ric{a2a2), where t] is the quantum efficiency 
of photodetection, and c is the velocity of light (we have put the charge of electron equal 
to unity so that the photocurrents are measured in number of electrons per second). For 
simplicity in what follows we shall consider the situation of rj = 1. Using Eq. ()3.10|) we can 
write the mean photocurrent (ii) measured in the transmission branch of the PBS as 

(ii) = {ii{ip, 6)) = -rjc So + 5*1 cos 2(^9 + (5*2 cos6' + 5*3 sinO) sm2ip , (3-11) 

where are the classical Stokes parameters. 

Equation (jH.llj) is the well-known formula for measuring the four classical Stokes parame- 
ters. The first three of them are obtained by removing the compensator [Q = 0) and rotating 
the transmission axis of the PBS to the angles (f = 0°, 45°, and 90°, respectively. The fourth 
parameter, S3, is measured by using a compensator with 6 = 90° or so-called quarter- wave 
plate, and setting the transmission axis of the PBS to (f = 45°. The four photocurrents are 
found to be, respectively, 

(^1(0°, 0°)) = ir^c(^o + ^i), 
(zi(45°,0°)) = lr/c(5o + 52), 
(2i(90°,0°)) = ^vciSo-S,), 
(^i(45°,90°)) = ^r]ciSo + S;). (3.12) 

Solving Eq. ()3.12|) for S*^ we can obtain all classical Stokes parameters from these four 
measurements. 



C. Observation of the fluctuation spectra of the quantum Stokes parameters 

In quantum optics in addition to the mean values of the quantum Stokes parameters 
(Sfj) their quantum fiuctuations are also taken into account. In this paper to describe 
the quantum fiuctuation we shall introduce the fiuctuation spectra of the quantum Stokes 
parameters. 
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Let us split the quantum Stokes operators S^(t) given by Eq. ()3.2|1 into the stationary- 
mean value = {Sfj} and small fluctuation 6S^{t), 

S,{t) = S^ + 5S^{t). (3.13) 

Taking the Fourier transform of 6S^(t), 

+00 

6S^{n) = ^ / SS^it)e''''dt, (3.14) 



/27r 

—00 

we can introduce the normally-ordered spectral correlation functions of the fluctuations 
6Sf^{Q) similar to the spectral correlation functions of the quadrature components in 
Eq. (Em, namely, 

(: 6S^{Q)6S,iQ') :) = {6S^6S,)n6{Q + Q'), (^ ^ u). (3.15) 

Here {6S'^)q are the spectral densities of the corresponding fluctuations and {6Sfj,6Su)n their 
cross-spectral densities. The symbol : . . . : means normal ordering of operators. 

To measure the spectral densities (55^)^ and the cross-spectral densities {6Sfj_6Siy)n of the 
quantum Stokes parameters given by Eq. ()3.15|) we can use an experimental setup similar to 
one that we have used for the measurement of the classical Stokes parameters (see Fig. 3). 
The difference is that instead of detecting the mean photocurrents (ii) and (^2) after the 
PBS, one observes now the photocurrent fluctuation spectra {6ip)^,p =1,2 defined as 

+00 

{6tl)n = J dt e*™((5v(0)5zp(t)), (3.16) 

—00 

where {6ip{0)6ip{t)) is the correlation function of the photocurrent fluctuations 6ip{t) = 
ip — (ip), and (ip) is the mean value of the photocurrent. Alternatively, one can add and 
subtract the individual photocurrents in the secondary channels and to investigate the sum 
i+{t) = ii{t) + i2{t) and the difference i-{t) = ii{t) — i2{t) of two photocurrents. In this 
case the information about the fluctuation spectra of the quantum Stokes parameters is 
contained in the fluctuation spectra 

+00 

{6tl)n= /rfte'™(5«±(0)5«±(t)). (3.17) 
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The photocurrent fluctuation spectra {Sip)Q and {Si'^)fi can be easily expressed through the 
spectral densities {SS'^)q and the cross-spectral densities {SSf^6S^)Q of the four quantum 
Stokes parameters. The results are conveniently presented in terms of the following linear 
combination of the three Stokes operators, S'l, S2, and S3, 



S = Si cos 2(f + {§2 cos 6 + S3 sin 6) sin 2if, 



(3.18) 



which is sometimes called a polarization observable [29|, |30|. We obtain the following ex- 
pressions for the fluctuation spectra {Sip)^ and {61"^) normalized to the shot-noise levels, 



1 + 



1 + 



2(ni) 



{SS^,)n + 2{SSoSS)n + {SS^ 



2(^2) 



(3.19) 
(3.20) 
(3.21) 
(3.22) 



where the corresponding spectral densities and cross-spectral densities of are defined ac- 
cording to Eq. (j3.15p . Here = {ii) + {12) is the shot-noise level of the photocurrent 
sum and difference, {rii) = (a|di), and (^2) = (0.20-2) are the mean photon numbers in the 
corresponding secondary channels after the PBS, and (n) = {rii) + (^2). 

Equations ()3.19|) - ()3.22|) are analogous of Eq. ()3.11|) for measuring the spectral densities 
of the quantum Stokes parameters. It is clear from these equations that with proper choice 
of angles 6 and if all nonzero spectral densities and cross-spectral densities of the Stokes 
operators can be measured. 



D. Relations between the spectral densities of the quantum Stokes parameters and 

of the quadrature components 

In Sec. II D we have provided analytical results for the fluctuations of the quadrature 
components (5X^(fi), 6Xy{Q), 6Yy{Q), and for their spectral densities and cross-spectral 
densities [see Esq. (I2.28p - (j2.33|) ]. Now we shall express the spectral densities of the quantum 
Stokes operators through the spectral densities of these quadrature components. As before, 
we shall restrict ourselves to the case of the x-polarized stationary solution when {rix) = 2Q^ 
and (uy) = 0. 
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Using the same normal rule of correspondence between the operators and their c-number 
representations as in Ref. Q we shall introduce the c-number variables S^{t) corresponding 
to the quantum Stokes operators S^(t). Since in Eq. ()H.2j) the Stokes operators are normally 
ordered, the same relation holds true for S'^(t) and the c-number variables ai(t) and a*(t), 
i = x,y. 

Linearizing the c-number variables S^{t) around their stationary values as 



S,{t) = S^ + 6S,{t), 



(3.23) 



we can express the fluctuations SS^lt) through the fluctuations of the field components 

6arc{t) and 6ay{t), 

6So{t) = 6Si{t) = V2Q[6a,{t)+6a:{t)), 
6S2{t) = V2Q[6ay{t) + 6a;{t)), 

SSsit) = -V2tQ[5ayit)-5a;it)). (3.24) 

Taking into account Eq. ()2.16|1 we obtain the following results relating the spectral densities 
of the Stokes operators with those of the quadrature components, 











(sslh = 




{5S25S3)n = 


8Q^i5Xy6Yy)n 



With help of these relations we arrive at. 



1 + 8k 
1 + 8k 
1 + 8k 



cos^ ^{6X^)n + sin^ ^{6X^)n 
sin^ ip{6X^)n + cos^ ^iSX^)n 



cos^ 2^{6X^)n + sin^ 2^{6X^ 



8)Q 



1 + 8k{SX^] 



To simplify Eqs. ()3.26p - ()3.28|) we have introduced the following shorthand notation, 

6Xein) = cose6Xy{n) - sine 5Yy{n), 
with its spectral density {6Xg)Q given by, 

{6X^)n = cos^ 9{5X^)n - 2 sin ^ cos e{6Xy5Yy)n + sin^ 9{5Yy%. 



(3.25) 

(3.26) 
(3.27) 
(3.28) 
(3.29) 

(3.30) 

(3.31) 
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The mean values of the individual photocurrents (ii) and {12), and of the photocurrent sum 
= {ii) + (^2) are equal to 

{ii) = 2g2/tcosV, (^2) = 2Q^Ksm^!f, = 2Q^k. (3.32) 

In the next section we shall investigate in detail the spectral densities of the quantum Stokes 
parameters and their cross-spectral densities. 



IV. POLARIZATION STATES OF LIGHT IN VCSELS 



A. Polarization squeezing 

The spectral densities {SS'^)q of the quantum Stokes parameters can be measured using 
any of three Eqs. ()3.19|) - ()3.2H) . Here we shall use Eq. ()3.21|) corresponding to observation 
of the noise spectrum {6i'^)Q{ip,9) of the photocurrent difference. With help of Eq. ()3.18|) 
we can bring the photocurrent noise spectrum {6i'^)Q{ip,d) to the form 

- {6 S2S S3) n 2 sin e cose +{6 Si) nsin^]}. (4.1) 

In this equation we have explicitly indicated the dependence of the observed noise spectrum 
on the angle 6 introduced by the compensator and the angle (p of the polarization beam 
splitter. 

The spectral densities {6Sq)q = {SSDq and (5S'|)f7 of the Stokes parameters Sq, Si and 
5*2 are measured by removing the compensator {6 = 0) and setting the transmission axis of 
the PBS to the angles (p = 0° and ip = 45°. The spectral density of the parameter 5*3 is 
obtained by using a compensator with 9 = 90° (quarter-wave plate), and setting ip = 45°. 
The corresponding photocurrent fluctuation spectra are given by, 

(5^)^.(0°, 0°)/(z+) = l + %{6S!)n, (4.2) 



{6^^)ni^5°, 0°)/(z+) = 1 + ^{SS^n, (4.3) 

2k 



{6zl)n{45\ 90°)/(h) = 1 + (4.4) 

In Fig. 4 we have shown the photocurrent fluctuation spectra given by Eqs. ()4.2|) - ()4.4|) for 
physical parameters close to that used in experiment namely, k = 100 GHz, 7 = 1 GHz, 
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7_L = 1000 GHzj'fs = 50 GHz, ujp = 40 GHz, a = —3, r = 6, and p = 1. The parameter 
describing the dichroism of the laser crystal, was set equal to zero in Fig. 4a, to = 10 GHz 
in Fig. 4b and to Ha = 50 GHz in Fig. 4c. 

Let us first discuss the case without dichroism (Fig. 4a). As seen from Fig. 4a, the 
spectral density {SSf)^ of the Stokes parameter 5*1 has a peak at a characteristic frequency 
Qi, while two other spectra (^S^j^ and {6S^)n for the Stokes parameters 5*2 and 5*3 exhibit 
peaks at another (higher) characteristic frequency fl2- These peaks are well-known from 
the theory of solid-state and semiconductor lasers and have their physical origin in the 
relaxation oscillations due to a periodic energy exchange between the active medium and 
the laser radiation. Since in our case there are two upper levels |a+) and |a— ) in the active 
laser medium, we have two subsystems where the periodic energy exchange takes place 
independently. First subsystem is described by the total population D of the upper levels 
and the Stokes parameter Si [see Eqs. ()2.17p ]. and its frequency of the relaxation oscillations 
is equal to Qi. In the second subsystem the relaxation oscillations take place between the 
population difference d and the two Stokes parameters 5*2 and 5*3 at the frequency ^2 [see 
Eqs. flTTR|l ]. 

Second important feature that one can observe in Fig. 4a is reduction of the quantum 
fluctuations of the Stokes parameter Si below the standard quantum limit at low frequencies 
Q in the case of regular pumping, p = 1. Thus, we can speak of phenomenon of polarization 
squeezing with respect to 5*1 in VCSELs with regular pumping. This result is to be expected. 
In fact, as follows from Eqs. ()3.2|1 . for the x-polarized stationary solution the Stokes parame- 
ter 5*1 coincides with the total number of photons in the laser field. It is well known from the 
literature 3] that a regularly pumped two-level laser can exhibit the sub-Poissonian photon 
statistics, i. e. the fluctuations of its photon number could be reduced below the standard 
quantum limit. One could therefore say that the polarization squeezing with the respect 
to Si in a regularly pumped VCSEL is the consequence of the sub-Poissonian statistics of 
photons. 

However, it is worth noting that the relation between the sub-Poissonian statistics of 
photons and the regular pumping statistics in VCSELs is not so direct as in the case of a 
two-level laser considered in j3|. Indeed, due to the degeneracy of the upper laser level on 
two sublevels \a+) and |a— ), the regular pumping of the total population D of the upper 
level remains random for each individual sublevel due to the partition noise. It turns out 
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that in the case of x-polarized stationary solution this partition noise does not contribute 
to the fluctuations of the total photon number and of the Stokes parameter 5*1. The reason 
for this is that, as follows from Eqs. ()2.17|1 . the fluctuations of the Stokes parameter 5*1 are 
coupled only with the fluctuations of the total population D and not with fluctuations of 
the populations of individual sublevels. 

The role of dichroism is illustrated in Fig. 4b and 4c. As seen from these figures, appear- 
ance of dichroism in the system has two major consequences. Firstly, the quantum noise 
reduction below the standard quantum limit in the spectral density {SSDq, of the first Stokes 
parameter is deteriorated by the factor + Kq). This deterioration has a clear physical 
explanation. Nonzero dichroism introduces random losses of the laser radiation inside the 
resonator at the rate Kq. The total decay rate of the laser field inside the resonator is now 
given by k + Kq, while the outcoupling rate determined by the transmission of the cavity 
mirror is equal to k. 

The second consequence of dichroism in the system is suppression of the relaxation os- 
cillations at the frequency ^2 related to the Stokes parameters 5*2 and 5*3. We can see from 
Fig. 4b that for small values of Ha {i^a = 10 GHz while k = 100 GHz) the peak of relaxation 
oscillations at ^2 becomes more pronounced. This is explained by the fact that for these 
values of we approach closer to the instability region. However, with increasing as in 
Fig. 4c the relaxation oscillations at ^2 rapidly disappear. 

The three spectral densities {6Sf)n, and {6Sf)Q in Fig. 4 can be also interpreted 

in terms of the uncertainty ellipsoid that we have mentioned in Sec. Ill A. Since the spectral 
densities depend on the frequency Q, one has to speak about the frequency- dependent un- 
certainty ellipsoid with tree major axis determined by the corresponding spectral densities. 
These spectral densities are normalized to the shot-noise level so that a sphere of unit radius 
in the Stokes-Poincare space corresponds to the standard quantum limit realized for a co- 
herent polarization state. As follows from Fig. 4a, for example, for a polarization-squeezed 
state in the area of low frequencies, where {SSDq is below the standard quantum limit, the 
uncertainty ellipsoid has the shape of a pancake. Instead, in the vicinity of the frequency 
of relaxation oscillations Qi this uncertainty ellipsoid takes a cigar-like shape with {6Sf)fi 
larger than two other components. 
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B. Cross-correlation spectrum of photocurrents 

Using the experimental setup shown in Fig. 3 one can also measure the cross-correlation 
function of fluctuations between the photocurrents ii{t) and i2{t), i. e. {Sii{0)6i2{t)) , or the 

corresponding cross-correlation spectrum of fluctuations, 

+00 

{Sii5i2)n = J dt e*™(5zi(0)5z2(t)). (4.5) 

—00 

Usually it is more customary to work with the normalized cross-correlation spectrum of the 
photocurrent fluctuations, 

Cuin) = , ■ (4.6) 

Using the Cauchy-Schwartz inequality one can demonstrate that this spectrum is normalized 
as |Ci2(^^)| < 1. Hence, Ci2(^^) = — 1 corresponds to the maximum anticorrelations between 
the two photocurrents, while Ci2(^^) = 1 to the maximum correlations. Experimentally this 
spectrum can be measured as. 

The normalized cross- correlation spectrum C^i^) can be expressed through the spectral 
densities and cross-spectral densities of the amplitude quadrature components 6X1 and 6X2 
as, 

Cuin) = , ^ 8.{6X,6X2h ^^^^^ ^4 



l + 8K{6X!)n^l + 8K{6Xl)n 
Using the relations between the field amplitudes di and 02 of the transmitted and reflected 

waves after the PBS and the incoming amplitudes and dy, given by Eq. ()3.1()|1 . we obtain 

{SXi5X2)n = cosipsmLp\{5X^)n~ {5X^)n\, 
{6Xl)n = cosV(<5^')c + sinV(5^e)n, 

{5Xi)n = sin2 ^i5X% + cos' ^{dX^h- (4.9) 

These relations allow us to express the cross-correlation spectrum Ci2(^^) in terms of the 
spectral densities {6X^)n and (5X|)n calculated earlier. 

In Fig. 5 we have plotted the cross-correlation spectrum Ci2{^) ioi (f = it /A and = 0. 
In this case the general result for Ci2{^) given by Eqs. ()4.8|) - ()4.9p is simplified to, 

Cx2{^) = ^ (4.10) 



4k 








1 + Ak 






n 
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Fig. 5a shows this cross-correlation spectrum for the case without dichroism and the same 
values of physical parameters as in Fig. 4. As follows from Fig. 5a, the cross-correlations 
are absent at high frequencies fl larger than 30 GHz. At lower frequencies of the order of 
15 GHz the curve of Ci2(^^) shows anticorrelations which turn to correlations at still lower 
frequencies of the order of 5 GHz. In the area of low frequencies Q smaller then 1 GHz one 
has again anticorrelations. 

This oscillating behavior of the cross- correlation spectrum Ci2{fl) is in full agreement with 
behavior of the fluctuation spectra of the Stokes parameters Si and 5*2 in Fig. 4a. Indeed, the 
cross-correlation function Cui^) is proportional to the difference of the spectral densities 
of the quadrature components {6X^)n — [or the corresponding Stokes parameters, 

i^SDn — (^S'Dq]. Therefore, for (5Xj)n > (5X^)n we have correlations between the two 
photocurrents, while in the opposite case - anticorrelations. 

Fig. 5b illustrates the same cross-correlation spectrum in presence of dichroism for differ- 
ent values of parameter Kq. As mentioned above, the essential role of dichroism is in the sup- 
pression of the relaxation oscillations. When Ka approaches the critical value Ka = 10 GHz 
of the instability border, the relaxation oscillations grow up and reinforce anticorrelations. 
Further increase of Ka results in suppression of the relaxation oscillations and respectively 
in transformation of anticorrelations into correlations for Ha larger than 50 GHz. 



C. Cross-correlations between the Stokes parameters 52 and S3 

For the x-polarized stationary solution that we consider in this paper, the linearized field 
operator E{t) from Eq. ()3.1|) can be approximately written as, 

= e^^*[y2Q + 6X,{t) + t6Y,{t)] [e, + -^{6Xy{t) + t6Yy{t))ey]. (4.11) 

This representation of the linearized field operator is very useful as it clarifies the physi- 
cal meaning of the quantum fluctuations of the four quadrature components that appear 
in Eq. ()4.1H1 . The fluctuations 6Xx{t) and SYxit) describe respectively the quantum fluc- 
tuations of the amplitude and the phase of the electromagnetic field -E(t). The quantum 
fluctuations of two other quadrature components 6Xy{t) and 6Yy{t) characterize the quan- 
tum fluctuations of the polarization of the field E{t). To see this more clear let us compare 
Eq. ()4.11|) with the classical expression often used in the literature on VCSELs (see for 
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example Ref. 0), 

^ e'^'\E\ [4 - (50 + i5x)ey] . (4.12) 

In this expression we have neglected the amplitude and the phase fluctuations of the fleld 
and have introduced the fluctuations 50 and 5x, 50 -C 1, 5x 1 of two angles and 
Xi that characterize the optical polarization state on the Poincare sphere. The flrst an- 
gle (0 < < tt) is called the polarization angle and it determines the direction of the 
polarization ellipse. The second angle x (^^/4 < X ^ ^/4) is the ellipticity angle. For 
x-polarized fleld both of these angles are zero. Comparing Eq. ()4.1H1 and Eq. ()4.12j) we 
conclude that these two classical fluctuations can be associated with their quantum counter- 
parts as 50 ——jJ^ and 5% — * — ^j^q - Taking into account Eq. ()3.24|1 we can also write 

'^'^2 , . 553 

50--^and5x--^. 

Thus, the quantum fluctuations of the Stokes parameter S2 characterize the fluctuations 
of the polarization angle, and those of the 5*3 - the fluctuations of the ellipticity angle. In the 
subsection A we have evaluated the fluctuation spectra of the Stokes parameters and S'3. 
However, as follows from Eq. ()3.25|) these two parameters are also cross-correlated. Hence, 
we shall introduce the cross-correlation spectrum C2^{^) between these two parameters in 
the same way as we did for characterization of the cross-correlations of two photocurrents. 

This cross-correlation spectrum is normalized as |C23(fi)| < 1 and can be experimentally 
determined from the measurements of the following three photocurrent fluctuation spectra, 

(5^)^,(45°, 0°)/(z+) = 1 + ^,iSSi)n, (4.14) 

(5^2_)n(45°,90°)/(z+) = l + ^(55|)n, (4.15) 

(5^)^(45°, 45°)/(z+) = l + -^[{6S^,)n + {6Sl)n + 2{SS2SSs)n]. (4.16) 

We have numerically evaluated the cross-correlation spectrum C23(fi) for the same values 
of physical parameters as in the previous subsection. In Fig. 6 we illustrate these spectra in 
the absence of dichroism (k^ = 0) and for two different values of equal to 10 GHz and 
50 GHz. 

As follows from this flgure, in the absence of dichroism the cross-correlation spectrum 
shows negative correlations at low frequencies Q less than 10 GHz. These anticorrelations 
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appear due to the coupling between the Stokes parameters 5*2 and 5*3 via the population 
difference d. For higher frequencies this coupling becomes less efficient and for VL higher 
than 30 GHz the ffuctuations of 5*2 and 5*3 become independent (C23 0). 

For nonzero dichroism the anticorrelations between S2 and 6*3 at low frequencies firstly 
disappear and then turn into positive correlations for larger values of Kq, for example at 

= 50 GHz. Thus, dichroism changes the nature of correlations between 5*2 and 5*3. 

V. CONCLUSIONS 

In conclusion we have presented a generalized and fully analytical theory of quantum 
fluctuations in VCSELs, proposed for the first time in Ref. Q]. The original results of our 
investigation are the analytical expressions for the spectral densities of the quadrature field 
components and of the corresponding quantum Stokes parameters. These analytical results 
facilitate the comparison between the theory and the experimental measurements. Moreover, 
we have included into the theory a nonzero linear dichroism of the semiconductor medium 
that was neglected in Ref. 

Our theory is very closely related to possible experimental observation of the quantum 
fluctuations in VCSEls that can be performed in a correlation-type measurement shown in 
Fig. 3. We have calculated analytically and illustrated graphically the typical fluctuation 
and cross-correlation spectra that could be observed in this type of measurements. Our 
theoretical results allow for direct comparison with experiments. 

We predict theoretically polarization squeezing in VCSELs when the quantum fluctua- 
tions of the Stokes parameter Si are reduced below the standard quantum limit. This phe- 
nomenon has its origin in regular pumping statistics of the active laser medium. However, 
the regularity in the pumping statistics alone is not sufficient for polarization squeezing 
in this type of lasers due to the partition noise between two upper sublevels in the laser 
medium. The second important feature of VCSELs that guarantees polarization squeezing 
is their dynamical behavior that couples the statistical properties of the Stokes parameter 
Si only with those of the total population of two upper sublevels. 

We have analyzed the role of linear dichroism and have concluded that it mainly influences 
the relaxation oscillations in VCSEls. These oscillations are typical for the solid-state and 
semiconductor lasers. The particularity of VCSEls is that in this case there are two types 
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of relaxation oscillations with clearly distinct characteristic frequencies fli and First 
oscillations (with frequency Qi) are related to the total population of two upper sub-levels 
and they contribute to the fluctuation spectrum of the Stokes parameter Si. The second type 
of relaxation oscillations (with frequency ^2) is connected with the population difference and 
its peak appears in the fluctuation spectra of the Stokes parameters S2 and 5*3. It turns out 
the dichroism dumps the relaxation oscillations of the second type and does not influence 
those of the first type. To understand this result let us recall that the relaxation oscillations 
appear in the lasers of the second type when the resonator losses are more rapid compared 
with those of the laser medium. As follows from Eqs. ()2.17j) and ()2.18|) dichroism increases 
the losses for the y-polarized light component coupled with the population difference d and 
does not change those of the x-polarized component related to population sum D. 
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FIG. 3: Experimental scheme for measurement of the spectral densities and cross-spectral densities 
of the quantum Stokes parameters. 
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FIG. 4: Photocurrent fluctuation spectra for the Stokes parameters 5*1, 5*2 and 5*3; a) without 
dichroism, Ka = 0, b) with dichroism, Ka = 10 GHz, and c) with Ka = 50 GHz. The values of 
other parameters are: k = 100 GHz, 7 = 1 GHz, j± = 1000 GHz, = 50 GHz, ujp = 40 GHz, 
a = 3 and p = 1. 
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FIG. 5: Cross-correlation spectrum Ci2{^) for ip = 7r/4 and = 0; a) without dichroism, = 
and b) with dichroism, Ka = 10 GHz and Ka = 50 GHz. The inset in a) illustrates the role of the 
statistical parameter p at low spectral frequencies. All other parameters are as in Fig. 4. 
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FIG. 6: Cross-correlation spectrum C23{^) without dichroism, = and with dichroism, 
10 GHz and Ka = 50 GHz for the same values of physical parameters as in Fig. 4. 



